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' — ' Abstract 
> 

Powerful new multi-kiloton liquid scintillator neutrino detectors, including NOi/A and LENA, 
will come on-line within the next decade. When these are coupled with a modest-power 
^ decay-at-rest (DAR) neutrino source at short-baseline, these detectors can decisively address 
ly-^ the recent ambiguous signals for neutrino oscillations at high Am^. These detectors are > 50 
O m long, and so with a DAR beam, the characteristic oscillation wave will be apparent over 
^ the length of the detector, providing a powerful verification of the oscillation phenomena. 
^ LENA can simultaneously perform i/^ — )■ appearance and z/g — )■ disappearance searches 
with unprecedented sensitivity. NOi/A is hkely hmited to z/g disappearance given its present 
design, but also has excellent sensitivity in the high Am^ region. For the appearance chan- 
ce nel, LENA could provide a stringent test of the LSND and MiniBooNE signal regions at 
> 5 cr with a reduced fiducial volume of 5 kt and a 10 kW neutrino source. In addition, 
the LENA and NOz/A disappearance sensitivities in z/g mode are complementary to the 
recent reactor anomaly indicating possible z/g disappearance and would cover this possible 
oscillation signal at the 3 a level. 
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1 Introduction 



Recent results from short-baseline (SBL) neutrino oscillation studies seem to point towards 
the existence of sterile neutrinos. The strongest indication comes from the LSND experi- 
ment (l}]4], which has reported a 3.8 o" excess of v^, events in a beam of v^. The associated 
Am^ of the assumed oscillation is too large to be explained with only three active neutrinos, 
and so oscillations involving sterile neutrinos are invoked which do not couple to W and 
Z bosons. The LSND result is supported by an apparent excess of z/g events in a beam of 

above 475 MeV observed in MiniBooNE [H] . This result is consistent with two- neutrino 

— ^ oscillations at 99.4% confidence level. 



More motivation has arisen from a recent reanalysis of the expected z/g flux emitted from 
nuclear reactors |6| that leads to an observed deficit of Ve at 98.6% CL. The reactor anti- 
neutrino flux prediction depends on accurate prediction of the reactor isotopes produced as 
a function of time convoluted with the spectra of anti-neutrinos produced by each isotope 
[T]. This anomaly arises from new calculations of the second ingredient, the anti- neutrino 
spectra, that update analyses from the 1980's [8 10 by including the latest information from 
nuclear databases and improving the calculational techniques [6]. The overall reduction 
in predicted flux compared to the existing data from SBL neutrino experiments can be 
interpreted as oscillations at baselines of order 10-100 m [IT] consistent with the LSND 
and MiniBooNE anti-neutrino results which require a mass squared difference of the order 
0.1-10 eV2. 

The SBL experimental results can be described using models which involve three active 



and one sterile (3+1) or two sterile (3+2) neutrinos 12 -21 , with relatively small mixing 
to the active flavors, allowing for oscillations with high Am^. While a sterile neutrino is 
theoretically well motivated with a large number of phenomenological consequences, the 
actual oscillation parameters required by the LSND and MiniBooNE anti-neutrino data are 
in conflict with various constraints imposed by other SBL neutrino oscillation experiments. 



most notably CDHS 22 , MiniBooNE disappearance measurements 23 and also the data 



from reactor experiments, like Bugey 24 



atmospheric and solar neutrino data in the required parameter range 18 



Also, there is no sign of sterile oscillations in 

This tension 



could be relieved by fitting neutrino and anti-neutrino results separately 25 
that neutrinos might oscillate differently than anti-neutrinos 



suggesting 



When analyzed within the 
context of the recent reactor anti-neutrino anomaly, a low, but acceptable, compatibility is 
found in the global fit using a (3+2) model 26 . The oscillation parameters found in these 
fits have Am^ ~ 1 eV^ associated with the sterile neutrino oscillations 

In the next five years, a series of small experiments will explore the (3+1) and (3+2) signal 
space further. This program includes continued MiniBooNE running in anti-neutrino mode 
and running of the MicroBooNE Experiment in neutrino mode 27 . Also, the upcoming 



KATRIN beta decay experiment is sensitive to a light sterile neutrino 28 . The Low Energy 



Neutrino Spectroscopy (LENS) detector, which is now in the prototype stage, may also 
run with radioactive neutrino sources in order to explore the question of sterile neutrinos 
29,30 . However, none of these experiments are expected to provide confirmation at the 



5 cr confidence level in the near future. 
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The next major step in the search for SBL oscillations can be achieved by pairing one of 
the ultra-large liquid scintillator detectors planned for the near future with a low energy 
neutrino beam. These experiments can study z/^ — )■ Ve and Ve oscillations. The high- 
sensitivity z/g disappearance search will bring unique information to the fits for comparison 
to the reactor z/g disappearance results. A relatively low power neutrino source producing 
neutrinos via the pion decay-at-rest (DAR) chain, is ideal. The DAR chain leads to a 
beam dominated by neutrinos between 20 and 52 MeV, with a well-defined flavor content 
of z^e, t'/i and z/^, as shown in Fig. [T] This source may be provided by a low energy proton 
accelerator with a beam impinging on a target/dump. Potentially, this can be a prototype 



for the cyclotrons planned for the DAEDALUS CP- violation search 31 33 . This would be 
a small, relatively inexpensive proton source which can be easily positioned within 20 m of 
the detectors. 

The energy range of the DAR beam is well suited to observe the L/E dependence of the 
oscillation wave across the length scales of presently planned detectors. A precise search for 
appearance and disappearance oscillations can be achieved by fitting the observed events 
with respect to expectation as a function of L/E. Besides observing oscillations from a 
simple counting analysis, these types of experiments also will confirm that oscillations are 
taking place by seeing the L/E variations within the detector. For a 40 MeV neutrino energy 
and Am? = 2 eV^, the characteristic oscillation length is Lose — 50 m which matches with 
the length of the liquid scintillator detectors being proposed for the near-future program. 
Because of the falloff of the neutrino flux with distance, the part of the detector 

within 50 m of the DAR source provides the main oscillation search sensitivity. Therefore, 
maintaining at least 50 m of detector length near the source appears to be best for covering 
the oscillation space of interest. 

In this paper, we consider two examples of scintillator detectors. Our example of an unseg- 



mented detector is LENA 34 , under consideration within the LAGUNA 35 project and 



the forthcoming LAGUNA-LBNO 36 design studies. The conclusions should be similar for 



other unsegmented detectors such as Hanohano 37,38 . Our example of a flnely segmented 



detector is NOz/A 39 , which is under construction for long-baseline studies using a beam 
from Fermilab. This detector is found to be less powerful than LENA, but has the advantage 
that it is already under construction and has an ideal space for the neutrino source. Each 
of these experiments can be used as designed, with only the addition of the beam source. 

We note that the large liquid scintillator detectors represent one of three types of future 
neutrino detectors under discussion. The other two types of detectors are water and liquid 



argon (LAr) based detectors. Ref. 40 provides an initial exploration of SBL oscillation 
searches using a DAR neutrino beam and the existing Super-Kamiokande detector with the 
addition of Gadolinium. The opportunities for an LAr detector will be explored in |4l] . 

This paper is organized as follows. We begin with the description of neutrino source and 
DAR flux in Sec. 2. In Sec. 3, we deal with the relevant SBL neutrino oscillation proba- 
bilities. After that, we describe the characteristics of liquid scintillator detectors (NOz/A 
and LENA) in detail in Sec. 4. We also discuss the possible charged current (CC) inter- 
actions of a DAR beam in liquid scintillator. In Sec. 5, we present our results for sterile 
oscillation searches using DAR-LENA and DAR-NOz/A setups. We summarize and draw 
our conclusions in Sec. 6. 
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Figure 1: The energy distribution of different flavor of neutrinos in a DAR beam. 

4 X 10^^ per year, per flavor (z/^, and Ue), 
1.6 X 10^® per year of z/g (4 x 10~^ compared to other flavors); 
Delivered as 100 kW average power, with 200 kW instantaneous power, 
(50% duty factor allowing equal beam-on and beam-off data sets); 
800 MeV protons on target; 
±25 cm smearing (assumed flat) on neutrino production point; 
20 m distance from average production point to face of detector fiducial region. 

Table 1: Characteristics of the neutrino source assumed in calculations, unless explicitly 
noted in the text. 



2 The Neutrino Source and DAR Flux 



In a stopped pion source a proton beam of ~ 1 GeV energy interacts in a low-A target 
producing 7r+ and, at a low level, tt" mesons. The pions then are brought to rest in a 
high-A beam stop. The 7t~ will capture. The 7r+ will produce the following cascade of 
decays 



TT 



1^ + 



1-^ e"^ + z/g + Uf, 

resulting in z/^, z/^ and z/g, but no Ue, as shown in Fig. [T] The resulting flux is isotropic. 

Certain choices of proton beam kinetic energy and neutrino beam stop design are crucial 
to producing the purest, highest rate DAR flux 



42 , 43 . As discussed in Ref. 33 , the 



production of DAR neutrinos is relatively flat and maximized in the region of beam kinetic 
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energy 700 to 1500 MeV for a given beam power. Selecting the lower energy range reduces 
the 7r~ production, which has a threshold of 602 MeV. Above 1500 MeV, kaon production is 
above threshold. It is for this reason that high energy proton beams, such as the 3 GeV beam 
proposed for Project X |44j, are poor sources for a DAR beam. The most efficient DAR 
production is on a light target which has tightly bound neutrons. In high-A targets, energy 



is wasted in neutron production 45 . High-A targets also have increased vr" production 
compared to light targets |46|. For this reason, a spallation neutron source, which uses 
high-A targets for neutron production, is not an ideal venue for the kind of experiments 
described here. In the optimized DAR beam designs, the light target for production is 
typically embedded in a high-A, dense material. This increases the probability that the 
7r~ produced in the light target will be captured before decay-in-flight (DIF), minimizing 
the decay chain that produces Ve backgrounds. Upstream targets, as are used in neutron 
production and also for isotope production, must not be used for this study, because these 
will substantially increase DIF backgrounds [sj. As a model of a DAR source, we use the 



DAE5ALUS design 33 



The DAE5ALUS accelerators are cyclotrons 47 -49 , an ideal low-cost source for low energy 



(800 MeV) protons. The one caveat to the use of a cyclotron as the accelerator is that the 
bunch spacing is typically a few tens of nanoseconds {e.g. one DAEDALUS design operates 
at 66 MHz [49]), hence much smaller than the muon lifetime. As a result these machines 
effectively operate as a continuous source (often called a "CW" source). It is straightforward 
to run this CW beam for milliseconds and then turn the beam off; however, implementing 
shorter spills is costly. DAEDALUS designs typically run the beam for ~100 ms periods. The 
required long spill precludes use of beam-timing to identify flavors, as has been suggested for 



DAR studies at spallation sources that employ expensive linacs and buncher rings 50 52 

The DAE5ALUS accelerators are being designed to produce more than 1 MW average power. 
The instantaneous power is a product of the energy and current of the protons on target, 
while the average power also accounts for the duty factor. The DAE5ALUS power needs are 
driven by the multi-kilometer distances required for the CP- violation search. The extracted 
beam is defocussed and the carbon target has a tapered entrance to spread the beam, to 
allow high instantaneous power on target. In order to cycle between sites and to allow for 
beam off running, the DAEDALUS experiment employs 10% to 20% duty factors. This leads 
to high instantaneous power requirements that are not needed in the SBL experiments we 
describe here. 

As will be shown in Sec. |5} an average beam power of 100 kW is sufficient to achieve most 
goals. We propose to allow some beam-off periods so that cosmogenic backgrounds can 
be measured and subtracted. A reasonable run-plan could have a 50% duty factor with 
a 200 kW instantaneous power on the target, delivered with the same time structure as 
planned for DAE5ALUS. Table [T] provides an overview of the assumptions used for the 
calculations in Sec. |5} unless otherwise noted. The low power requirement, with similar 
time structure to DAEDALUS, would make these accelerators ideal as prototypes for the 
high-power DAEDALUS machines. 

The position of the neutrino production point depends on the proton interaction length in 
the target material, the pion stopping length in the target, and shielding and tapering of 
the target entrance to spread the beam power. We will assume a ±25 cm smearing for 
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the neutrino production point in this study. We take this to be flat in L, which is an 
overestimate of the smearing. We also assume that, on average, neutrinos are produced 20 
m from the detector face to allow for accelerator shielding. We assume DAE(5ALUS-level 
contamination, which is at 4 x 10~^ of the z/g in the beam. 

While the ratio of flavors and the energy dependence of the flux are well understood in a 
DAR flux, the overall normalization is not well predicted. The dominant error comes mainly 
from the uncertainty on the pion production per proton on target 42 , 43 , 46 . We assume a 



10% correlated normalization error on all flavors, which is slightly more conservative than 
We also assume a 20% error on the tt" DIF background [3j. 



m 



42,43 46 



3 (3+n) sterile neutrino oscillation hypotheses 

If one assumes CPT invariance and no matter effects, then the probability for a neutrino 
produced with flavor a and energy E, to be detected as a neutrino of flavor (3 after traveling 



a distance L, is 53 ,54 



P{i^a ^ i^p) = - 4 Ei>, n{U*^^UpiU^,Ul^) sin^ Xi, + 

2 E.>, X{Ul,Up.U^,U;^) sin 2x,,- (1) 

where TZ and X denote the real and imaginary parts of the product of mixing matrix elements, 
respectively. In Eq. [T| = e,/x,r, or s, {s being the sterile flavor); z, j = 1, . . . ,3 + n 
{n being the number of sterile neutrino species). The neutrino mass splitting is given by 
Am^j = ml — in eV^. The L/E dependence associated with oscillations appears within 
the term Xij = AmfjL/4E, where L is in m and E is in MeV. The U represent the elements 
of the mixing matrix that connects the mass to the flavor eigenstates. For anti-neutrinos, the 
oscillation probability follows Eq. [TJwith the replacement of U with its complex-conjugate 
matrix. Therefore, if the elements of the mixing matrix are not real, neutrino and anti- 
neutrino oscillation probabilities are not identical. Imaginary parameters enter the mixing 
matrix through CP-phases, which, in turn, lead to differences in the neutrino versus anti- 
neutrino oscillation probability. 

In a (3-1-2) model, two sterile neutrinos are added at the eV scale with the three active 
neutrinos. In the SBL approximation where the largest mass-splittings dominate and 
Amli ~ ^"^31 ~ 0, following Eq. [l| the z/^ Ve appearance oscillation probability is 
given by 

P{y^^Vf>) = 4|[/e4nt/^4psin^X4i 
+ 4|f/e5p|f/^5psin^X5i 

+ 8\Ue4Uf,4Ue5U^5\smX4lSinX5iCOs{x54 + 5) (2) 

where 6 = argiUlJJf^JJebU*^) is the CP-phase. We use this expression to estimate the signal 
event rates for appearance studies (see Sec. 5) in (3-1-2) models for the best fit parameter 
values as given in Table |2} We can see that the probability in Eq. |2] depends on two 



5 











\Ue,\ 










(5/7r 


A 


Ref. 


26 


0.47 


0.128 


0.165 


0.87 


0.138 


0.148 


1.64 


B 


Ref. 


25 


0.39 


0.40 


0.20 


1.10 


0.21 


0.14 


1.1 



Table 2: The 1st row depicts the parameter values at the global best fit points for 3- 
model as described in 



26 



25 which has 



The 2nd row shows the 3+2 best fit values from 
been derived using all the data sets of SBL appearance experiments. Here mass splittings 
are shown in eV^. 



independent mass splittings {Aml^, Aml^) and on the combinations \Ue4,U^4^\ and If/est^^isl- 
Therefore, including S, the total number of independent parameters in this channel is 5. 

In a (3+1) model with only one sterile neutrino, Eq. [2] takes the form 

P{i^^l i^e) = 4:\Uei\^\Uf,i\'^ sin^ = sin^ 29^e sin^ 2:41 (3) 

where sin^ 26'^e is the effective mixing angle. One can see that CP-phase does not appear 
in the (3+1) case. 

Again following Eq. [T] with the SBL approximation, the — )■ I'e disappearance oscillation 
probability in a (3+2) model can be written as 

PiUe^Ue) = 1 -4(1 - |t/e4|^ - |f/e5n(|f/e4rsin2x41 + l^/esTsin^Xsi) 

- 4|f/e4nf/e5psin2x54, (4) 

which we use to simulate the survived event rates in (3+2) models. In this channel, we have 
4 independent parameters. 

In a (3+1) model, Eq. |4] simplifies to 

P(z/e — )■ Z/g) = 1 — 4|[/e4| (1 " \Ue4\ ) siu X41 = 1 — siu 26'ee siu 0:41, (5) 

where sin^ 29ee is the defining mixing angle for SBL electron neutrino disappearance. 



4 Liquid Scintillator Detectors 



Scintillating-oil detectors use well understood technologies, but the planned next steps in- 
volve significant increases in scale. Segmented scintillator detectors have been traditionally 
used for GeV-scale experiments. The largest pure-scintillator, finely-segmented example is 
the 170t ton BNL 734 experiment 



mo 



57 



55 . On the other hand, KamLAND 56 



and Borex- 

have demonstrated ton-scale unsegmented liquid scintillator detectors, with the 
latter demonstrating a very low level of background, allowing low-energy (~ 1 MeV) neu- 
trino studies with low systematics. 

The next generation of these detectors are more than an order of magnitude larger than 
those described above. NOz/A, which is under construction at Ash River, Minnesota, US, 
is a segmented detector. This detector will come on-line in 2013 [58]. LENA, which is an 
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Figure 2: Cross-sections of all the flavor selective processes used in this study. 

unsegmented, very-low background detector, is proposed for tlie LAGUNA |35j project in 
Europe. Because site-selection for LAGUNA is only just underway, this detector will come 
online later than NOz/A, around 2020. 

Both NOz/A and LENA are intended for long baseline oscillation studies. The SBL running 
does not conflict with the long baseline running for several reasons. First, the long baseline 
neutrino beam has energies > 100 MeV, producing events that are easy to separate from 
the lower energy DAR events. Second, the long baseline beam timing allows one to gate out 
the DAR beam for a few milliseconds around the the long-baseline few-microsecond spill. 



4.1 DAR beam interactions in liquid scintillator oil 

As discussed above, a DAR beam consists of i^^, z/^, z/g with a small admixture of z/g with 
energies ranging up to 52.8 MeV. Because of the low beam energy, a muon cannot be 
produced in a CC interaction. Therefore, detectors with oil-based targets (CH„) are limited 
to observe CC interactions involving the electron flavor scattering. Inverse /3-decay, (IBD), 



z>e + P — ^ e"*" + n 59 is the golden channel for these detectors to detect the appearance of 
Vf.. This process has a very low kinematic threshold of 1.81 MeV and provides a very useful 
delayed coincidence tag between the prompt positron and the delayed neutron capture by 
a proton, n + p — )■ d + 7 (2.2 MeV) after an average time of 250 /xs. 



Another CC reaction that we consider is z/g+^^C — )■ e +^^Ng.s. 60 with a relatively high 



kinematic threshold of 17.33 MeV. This process is used to perform the disappearance 
searches using z/g. In this cross-section, we only consider the contribution coming from 
the transition to the ^^N ground state, which can be identified by the detection of the 
prompt electron, followed within a 60 ms window by the positron from the beta decay of 
the ^^Np^s with a mean lifetime of 15.9 ms. The cross section for this exclusive reaction 
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Detector 

V^Ildl dCtei lotlCo 


NOz/A Far 
jjetectoi 


LENA 


Shape 


Rectangular 


Cylindrical 


Fiducial Mass 


14 kt 


(5 - 50) kt 


Overburden 


3 m earth-equivalent 


1450 m of rock/4060 mwe 
(Q3 ryhasalmi 


Solvent 


CH2 


LAB (C18H30) 


Neutrino Energy 
Threshold 


38 MeV (Dis) 


20 MeV (App) 
33 MeV (Dis) 


Detection Efficiency 


50% (Dis) 


90% (App) 
80% (Dis) 


Energy Resolution, cr(E) 
in units of MeV 






15%v/E/MeV 61 


10%v/E/MeV 62 


Signal error (syst.) 


20% (Dis) 


10% (App) 
15% (Dis) 


Background error (syst.) 


5% [Non-beam] (Dis) 


20% [Intrinsic z/g, 4 x 10"^] (App) 



Table 3: Detector characteristics used in the simulations. "App" applies to the appearance 
analysis. "Dis" applies to the disappearance analysis. 



^^C(z/e, e )^^Ng.s. is well measured compared to the inclusive cross-section for transitions to 



excited states of ^^N and has an error of 5 to 10% 60 . Probably the contribution from tran- 
sitions to excited states of ^^N can also be used and would increase the event rate by about 
30 to 40%. The excited state cross-sections are less well known but since the power of the 
measurements presented here mainly comes from an L/E shape analysis, these uncertainties 
are not very important. 

Fig. [2] shows the cross sections both of the CC flavor selective processes used in this study. 
The cross section for i/e+^^C — )■ e"'"+^^B compared to that for IBD scattering is sufficiently 
low to be neglected from these studies. We also do not consider the neutrino-electron 
scattering interaction because the rates are two orders of magnitude lower than the rate of 
^^C scattering. 



4.2 A Finely-Segmented Detector: NOz/A 

NOi/A [39] is a segmented, scintillator detector constructed for long baseline neutrino os- 
cillation studies around 2 GeV, which is 40 times higher in energy than the DAR events. 
Because this detector is not designed for ~ 50 MeV events, it is not ideal for the studies we 
present here. Nevertheless, we show that NOi/A has the capability to produce significant 
disappearance results with a few years of DAR beam running. 

NOi/A consists of PVC cells filled with scintillator oil, assembled into a 15.7 x 15.7 x 67 
m'^ detector. The scintillator is mineral oil base doped with 4.1% pseudocumene [1,2,4- 
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Trimethybenzene] . The PVC structure (CH2— CHCl), which is dead region, represents 30% 
of the 14 kt fiducial mass. For simphcity, we assume that the target as entirely CH2 for 
this analysis. Individual cells are 3.9 cm high x 6.0 cm along the beam direction x 15.5 m 
long transversely and are glued in X and Y planes. The orientation of the planes alternates 
throughout the detector. 

Light produced in the scintillator is absorbed and re-emitted with wavelength shifting (to 
550 nm) optical fibers. These run the length of the PVC cell and double back to readout 
APDs. The APDs have a rather high noise rate but will be cooled to reduce the noise to be 
below the cosmic-ray muon rate. The attenuation length of the fiber is roughly the same 
as the 15 m length of the cells. Triggering requires at least one 1 hit in adjacent X and Y 
planes. The detector is below ground level with overburden shielding of 3m of earth. This 



modest shielding leads to a 20 kHz rate of through-going and 1 kHz of stopping muons 61 
in the detector. 

The NOz/A setup is ideal for the proposed cyclotron-dump design. The cyclotron can be 
located on the surface, allowing easy maintenance, while the target/dump can be located 
in an open area of the detector hall which is used for staging and not needed after the 
detector is constructed. The incoming beam will impinge on a target embedded in the 
dump. Targeting the beam perpendicular to the detector direction is an extra precaution 
against DIF backgrounds. In addition, the ~ 50 m length of the detector is appropriate for 
neutrino oscillation measurements with 20 to 50 MeV neutrinos probing the Am^ ~ 1 eV^ 
region. 

We expect that the shielding around the dump will lead to a very low fast neutron back- 
ground. Because the detector staging area is quite large, it will be possible to add up to 
10 m of extra shielding, if required. Beyond this, the detector is self-shielding. Therefore 
we do not expect a significant fast neutron background from the cyclotrons. It is for this 
reason that we can place the face of the detector at 20 m from the center of the target. A 
summary of the assumptions we make about the NOi/A detector in the analyses is given in 
Table E 

The fine segmentation of the detector allows determination of the vertex to a few cm. As a 
result, the smearing in the measured L value will be dominated by the uncertainty in the 
neutrino production point for which we use a ±25 cm fiat distribution. 

The disappearance study makes use of u^+^'^C — )■ e~+^^Ng.s. interaction with a visible 
reconstruction energy threshold of 20 MeV which corresponds to a neutrino energy threshold 
of 38 MeV. This threshold significantly reduces the accepted fiux for the measurements, 
as seen in Fig. [l] but still gives about 34k events per year (see Sec. [5]). We assume 50% 
efficiency for this process. 

The z/e scatters will be reconstructed along the 67 m length of the detector, allowing the 
oscillation wave to be fit as a function of L/E. Since the target nucleus is much heavier 
than the outgoing electron, the incident neutrino energy is simply related to the the visible 
outgoing electron energy independent of the scattering angle. The electron energy resolution 
is listed in Table [H 

The largest background source will be from the 10^° Michel electrons/year produced by 
stopped muon decay. These are produced with nearly the same energy dependence as the 
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electrons from CC v^, scatters. The Michel electron events can be identified and vetoed 
by tracking the parent muon. The entire target can be used to identify muons before a 
candidate z/g scattering event, which may require that a substantial part of the upper region 
of the target be used for this veto. In the analysis presented in the following section, we 
consider a range of veto capabilities giving from 10,000 to 50,000 total un- vetoed Michel 
background events per year. We also show the no-background case for reference. We assume 
that the 50% beam-off running can measure this background with a 5% error in each bin. 

In principle, the NOz/A detector can also tag z/g + p — e^-|-n events through the coincidence 
signal of the positron light followed by light from the 2.2 MeV photon emitted in the neutron 
capture on H. In practice, observation of the the neutron-capture signal in NOz/A will be 
quite difficult. First, the signal from the 2.2 MeV 7 is diffuse, being produced by electrons 
from multiple Compton scatters. Because the NOz/A detector has 30% dead material, there 
is a high probability that a substantial fraction of the 2.2 MeV energy deposit will be 
unobserved. Second, the attenuation length of the wavelength shifting fibers is on the scale 
of the 15 m length of the scintillator tubes. Thus, the detector is highly inefficient for very 
low light deposits. Third, the false coincidences with comic muon produced background may 
be very high in a surface detector. These are complex analysis issues that would require 
further study to quantify and we therefore choose not to present the appearance capability 
of NOz/A in this paper. 



4.3 An Unsegmented Detector: LENA 

The LENA (Low Energy Neutrino Astronomy) detector [34] is optimized for low energy neu- 
trino detection. The experiment is specifically designed for state-of-the-art measurements 
of supernova neutrinos (bursts and diffuse), which are in the same energy range as DAR 
flux, and solar neutrinos, geoneutrinos, and reactor neutrinos, all of which are an order of 
magnitude lower in energy. As a result, this detector is substantially more powerful than 
NOz/A for the DAR-based SBL program we consider here. 

The low energy neutrino studies require substantial shielding of at least 4000 mwe. The 
search for a possible underground site for LENA has been performed in the framework of 
the LAGUNA design study [35]. Two possible choices are Pyhasalmi and Frejus. Both the 
sites have rock shielding > 4000 mwe resulting a very low rates of through-going muons. 



at the level of ~ 5 x 10~^ Hz/m^ 34 . The stopping muon background rate is negligible. 
The rate of atmospheric electron neutrino interactions in the DAR energy range will also be 
negligible. Nevertheless, we assume that there will be 50% beam-off running. This allows 
cross-checks of the backgrounds as well as allowing the simultaneous data-taking for the 
supernova and SBL studies. 

The LENA detector is planned to be 100 m in length and 30 m in diameter. If a flducial 
volume of 13.6 m radius is chosen, then the detector will contain 50 kt of liquid scintillator 
while the outside region will be fllled with water to act as a veto for muons and shield for 
neutrons. The detector will have 30% coverage from a combination of direct photocathode 
coverage and light concentrators. 

In the studies below, we also consider several smaller flducial volumes. None are in the 
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present LENA plans. However, a 5 to 10 kt prototype may be attractive along the path to 
LENA and pressure at the base of the detector on the PMTs may lead preferance for a 25 
kt design. Therefore, we also include these possibilities in our discussion. 

Both CC ^^C events for disappearance and IBD scattering events for appearance can be 
used in the analysis. For LENA, we consider LAB solvent (CigHso) with 6.6 x 10^^/m^ free 
protons and 4.0 x 10^^/m'^ carbon nuclei in a 50 kt detector. Demonstration that oil can be 



purified well below the levels required for the SBL physics come from the KamLAND 56 
and Borexino [57] experiments. The specific assumptions about the LENA detector relevant 
to this study appear in Table [3j 

The arrangement of the cyclotron and target/dump will be designed after the location of 
the detector is established. Tunnel access where the neutrino source can be moved into the 
detector region is most attractive. We assume that the neutrino source will be located on 
the long-axis of the cylinder and will be 20 m away from the cylinder surface. 

LENA can explore both appearance and disappearance oscillations. The signals are the 
same as discussed for NOi/A: z/e + p — > e"^ + n and z/g+^^C — ?■ e~+^^Ng s., respectively. Be- 
cause the detector is designed for excellent reconstruction of geoneutrinos and astrophysical 
signals at 2 MeV, there is high efficiency for reconstructing the 2.2 MeV neutron capture 
7. We consider 90% detection efficiency for appearance studies and 80% efficiency for dis- 
appearance, although this is likely to be very conservative. Vertex reconstruction can be 



expected to have a 5 cm uncertainty 34 . This is negligible on the scale of the ±25 cm 
uncertainty from the neutrino source. 

We assume that there is very little background in LENA. The > 4000 mwe shielding reduces 
the cosmic ray background to a negligible level. The next largest natural background is from 
atmospheric neutrinos. These are at the <1% level for the 100 kW source. It should be noted 
that the atmospheric muon neutrino background is much lower than in a water Cerenkov 
detector. This is because, in a scintillation detector, all muons from charged current events 
are identified, whereas in the water detector some are below Cerenkov threshold. We assume 
a > 10 m wall of undisturbed rock (or equivalent steel or concrete shielding) in the design 
that will reduce beam backgrounds from neutrons produced at the accelerator to a negligible 
level at the detector. The additional water shielding further protects the detector. 

Disappearance searches proceed in the same manner as NOz/A. The signal is the variation 
of z/g+^^C —7- e^+^^Ng.s. across the length of the detector. For the oscillation sensitivity 
estimates, the simulated data, taking into account the detector energy resolution and the 
neutrino source position smearing, are binned according to energy and position in 65 equally 
sized L/E bins. A visible energy cut of 16 MeV is imposed for the events, which corresponds 
to a 33 MeV cut on the incoming neutrino energy. This visible energy cut will practically 
eliminate backgrounds from atmospheric and reactor neutrinos as well as environmental 
radioactivity. Normalization is included but because these errors are large, the fits are 
dominated by the L/E shape dependence. 

LENA is designed also to have high efficiency for z/g + p — )• + n (IBD) events at neutrino 
energies in the 10 to 50 MeV range identified by a delayed coincidence between the outgoing 
positron and capture of the neutron. This is the signal for diffuse supernova neutrinos, and 
an important process for supernova burst events - both of which are key physics goals of 
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(A : Ref. |26]) 


(B : Ref. [25] ) 


InTT'TFlQir' 7/ 

iiiiiiiibiL y ^ 
Background 


50 kt 


13.58 m 


100 111 


12985 


32646 


1450 


25 kt 


10.78 m 


79.37 111 


7787 


18356 


875 


10 kt 


7.94 111 


58.48 111 


3753 


7964 


443 


5 kt 


6.3 m 


46.42 111 


2080 


4044 


261 



Table 4: The expected number of signal and intrinsic beam background events in 5 to 
50 kt LENA detector. While varying the fiducial mass of the detector, we have kept the 
density and the aspect ratio (length/diameter) same in all the cases. In column 4 and 5, 
the signal events have been computed using the two different sets of parameter values in 
(3+2) model as given in Table [2] Here we have used total 4 x 10^^ z/^. The intrinsic z/g beam 
contamination is 4 x 10~^. 



LENA. Thus, LENA is specifically designed for lower energies in contrast to NOi/A. Using 
the IBD events, LENA can make a very precise search for z/g appearance by again binning 
and fitting events as a function reconstructed LjE. For our studies, the visible energy 
threshold is set to 19 MeV corresponding to a neutrino energy cut of 20 MeV. Again, this 
requirement will reduce the non-beam backgrounds to negligible level. 

This puts LENA in a unique position to measure both appearance and disappearance in 
the same detector through the characteristic change in the oscillation wave with L. That 
would be a very powerful signal in support of sterile neutrino models and offers the best 
opportunity to disentangle (3+1) from (3+2) oscillations. 



5 Results 

In this section, we first discuss SBL i/^ z/g appearance with LENA. This is key to confirm- 
ing or refuting the LSND (l}]4j and MiniBooNE [5] anti-neutrino results. Then we focus on 
SBL Vf. disappearance studies with LENA and NOz/A, which are complementary searches 
to the reactor anomaly z/g disappearance signal (6|[lT]. 



5.1 Appearance Mode 

The LENA appearance search uses the IBD signal. Following Ref. [34], we consider a cylin- 
drical unsegmented liquid-scintillator detector of 50 kt fiducial mass with 100 m in length 
and 13.58 m in radius as a reference choice. We also study the impact of smaller-sized LENA 
type detectors with fiducial masses of 25, 10 and 5 kt. In all cases, we keep the fractional 
pliotodetector coverage and aspect ratio (length/diameter) same as we consider for 50 kt. 
In addition, while varying the fiducial mass of the detector, we keep the detector charac- 
teristics unchanged as described in Table [3j The main source of background is the intrinsic 
Vf. beam contamination from vr" decays (see Table [g, because cosmogenic backgrounds are 



very low (see Sec. 4.3). The neutrino energy threshold for this analysis is 20 MeV which 
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Figure 3: The ratio of signal events estimated for a given set of oscillation parameter values and the 
events which have been computed assuming a hypothetical case where the v^^ — > v^^ oscillation probability is 
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The 



one is shown as a function of reconstructed LjE in 50 kt LENA. Here we have used total 4 x 10 
solid red line is for the parameter values in (3+2) model as given in the upper row of Table [2] The blue 
dashed line is computed for (3+1) best fit values: Am|j^ = 0.57 and sin^ 26*^2 — 0.0097 as given in j25j. See 
text for more details. 



renders potential backgrounds coming from supernova and radioactive decay negligible. We 
present all our results assuming a fiat resolution on the neutrino length distribution of ±25 
cm dominated by the uncertainty in the neutrino production point. 

In Table |4] we present the signal and intrinsic Ve background event rates in 5 to 50 kt LENA 
type detector including the source and the detector characteristics as shown in Table [T] and |3] 
respectively. In column 4 and 5, the signal events have been estimated for two different sets 
of parameter values in (3+2) model as given in Table [2] Here we consider a total neutrino 
flux of 4 X 10^^ z/^ from DAR beam. One should note that the number of signal events 
does not increase linearly with the fiducial mass because of the geometry of the detector 
and the isotropic nature of the flux. While increasing the fiducial mass from 5 to 50 kt, the 
length of the detector increases by a factor of 2.15 and since the source is located on the 
axis of the cylinder and 20 m away from the beginning of the detector, a very small amount 
of neutrinos reach the ends of the detector due to the 1/L^ suppression from the isotropic 
source. 

With a source-to-detector-face distance of 20 m, the accessible L range in 50 kt detector is 
20 - 120 m and with an energy range of 20 - 52.8 MeV in DAR beam, the L/ii^-dependence of 
the oscillation pattern can be well observed inside the detector. This is illustrated in Fig. [3} 
where in each reconstructed LjE bin, we show the ratio of signal events for a given set of 
oscillation parameter values to the number of events computed assuming a hypothetical case 
where the z/^ — )■ z/g oscillation probability is one. This has been done just to demonstrate the 
oscillation pattern inside the detector as a function of reconstructed LjE for a particular 
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r lU-X. 


r LU.X. 




(A : Ref. |26l) 


(B : Ref. [25] ) 


50 kt 


0.912 X 10^'^ 


0.302 X lO^y 


25 kt 


1.535 X 10^9 


0.539 X 10^9 


10 kt 


3.235 X 10^9 


1.27 X 10^9 


5 kt 


5.935 X 10^9 


2.6 X 10^9 



Table 5: It shows the amount of neutrino flux needed to exclude the two different sets 
of oscillation parameter values as given in Table |2] at 5 cr CL with 5 degrees of freedom 
{Ax^ = 37.09) using 5 to 50 kt LENA in appearance mode. 



set of parameter values. The red solid line corresponds to the parameter values in the (3+2) 
model as given in the upper row of Table |2| denoted by A. The blue dashed line has been 
drawn for the (3+1) best fit values {Aml^ = 0.57 eV^ and sin^ 29^e = 0.0097) obtained using 



LSND, MiniBooNE anti-neutrino and KARMEN 63 data sets in [25]. As one can see from 
Fig. [3} the demonstration of the oscillation wave is dramatic in the long LENA detector and 
can provide a powerful handle to discriminate between (3+1) and (3+2) schemes. 

For our statistical analysis, we bin our signal and background events into 94 equally sized 
L/E bins and consider uncorrelated normalization errors on the signal and intrinsic Ue 
background of 10% and 20% respectively. These uncertainties are fully correlated between 
the L/E bins and are included in the analysis using the pull-term method as described 



in e.g. Ref. 64,65 . For the fitting, we perform the usual x analysis using a Poissonian 



likelihood function. 

In Table [5] we present the amount of flux needed for a 5 to 50 kt LENA type detector to 
exclude the two different sets of oscillation parameter values in the (3+2) scheme as given 
in Table [2] The numbers presented here for 5 a CL (5 dof) under the assumption that we 
have only intrinsic z/g beam background. We can immediately see that with a very modest 
flux and a small size LENA detector, we can check these test points at high significance. 
This also shows that the L/E dependence is very important for rejecting the background 
and therefore reducing the sensitivity to systematic errors. The ability to observe the L/E 
dependence is crucial if a signal is observed, since it will allow one to establish or refute 
oscillations as the underlying physics explanation. 

Fig. |4] presents the sensitivity limit of the DAR-LENA setup to sterile neutrinos in the (3+1) 
model at 5(7 confidence level (2 dof) using the appearance mode. We compare our results 
with the allowed region at 99% CL (2 dof) from a combined analysis of the LSND and 
MiniBooNE anti-neutrino signals [26]. The results are presented for four different choices 
of the fiducial mass of the detector. The left panel shows the sensitivity for 4 x 10^*^ z/^ 
which can be achieved in one year with a small 10 kW average power machine. The right 
panel exhibits the same for our reference choice of 100 kW machine. It can be seen from the 
left panel that a 5 kt LENA type detector and a one year run at 10 kW average power is 
more than sufficient to exclude the parameter space suggested by the combined fit of LSND 
and MiniBooNE anti-neutrino data at 5 a CL. Note, that the sensitivity is limited by the 
magnitude of the beam background and therefore does not improve linearly with the size of 
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sin 20j,e sin 20^e 

Figure 4: Sensitivity limit of DAR-LENA setup to sterile neutrino oscillation in the (3+1) model at 5 cr 
CL (2 dof) using appearance mode. The green/gray shaded area is the allowed region at 99% CL (2 dof) 
from a combined analysis of the LSND and MiniBooNE anti- neutrino signals [26] . Left (right) panel shows 
the results for 10 (100) kW average power machine which can deliver 4 x 10^° (4 x lO'^^) in one year. In 
both the panels, results are shown for four different choices of the fiducial mass of the detector. 

the detector. In the right panel, the 50 kt detector has a reach up to sin^ 26 = 0.0001 at 
Am|i = 2 eV^ with 4 x lO^^ iy^. 

5.2 Disappearance Mode 

Both LENA and NOz/A offer the possibihty to study oscillations to sterile neutrinos in the 
disappearance mode by using the CC reactions of on ^^C. In Table |6] we present the total 
number of CC Ug scattering events on ^^C in 5 to 50 kt LENA type detectors using the 
information on the source and the detector characteristics from Table [T] and [3] respectively. 
Column 4 and 5 show the number of survived events after oscillation using the two different 
sets of (3+2) parameter values given in Table |2] The last column shows the total event rate 
without any oscillation. Here we have used total flux of 4 x 10^^ z/g with a neutrino energy 
threshold of 33 MeV and an efficiency of 80%. Table [T] shows the same for the 14 kt NOz/A 
far detector but with a neutrino energy threshold of 38 MeV and an efficiency of 50%. We 
can expect ~ thirty-four thousand events in NOz/A with no oscillation considering only the 
contribution of the transition of ^^C to the ^^N ground state [60]. From Table [6] we can 
see that the impact of disappearance on the total event rate is 6.3%/23.4% for the A/B 
parameter sets in the 50 kt LENA detector. One should note that this is not a mere counting 
experiment and that the L/E pattern of the oscillations provides most of the measurement 
sensitivity. This can be seen in Fig. |5] where we plot the ratio of events with and without 
sterile oscillation as a function of the reconstructed L/E m the 50 kt LENA detector. With 
an L range of 20 - 120 m and with an energy range of 33 - 52.8 MeV, the L/£'-dependence 
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50 kt 


13.58 m 


100 m 


170191 


139119 


181672 


25 kt 


10.78 m 


79.37 m 


102726 


85271 


109590 


10 kt 


7.94 m 


58.48 m 


52105 


43940 


55439 


5 kt 


6.3 m 


46.42 m 


30874 


26321 


32735 



Table 6: The total number of CC v^, scattering events on ^^C in 5 to 50 kt LENA. In 
column 4 and 5, the survived number of events have been computed using the two different 
sets of oscillation parameter values in (3+2) model as given in Table [2] The last column 
shows the total event rate without any oscillation. Here we have used total 4 x 10^^ z/g with 
an energy threshold of 33 MeV and an efficiency of 80%. 
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(B : Ref. (25)) 


Evts, No Osc 


14 kt 


67 m 


15.7 m 


15.7 m 


32388 


27407 


34415 



Table 7: The expected number of CC v^. scattering events on ^^C in 14 kt NOz/A far 
detector. In column 5 and 6, the survived events have been computed using the two different 
sets of oscillation parameter values in (3+2) scheme as given in Table [2] In the last column, 
we have the total event rate without any oscillation. Here we have used total 4 x 10^^ v^. 
with an energy threshold of 38 MeV and an efficiency of 50%. 



of the oscillation pattern can be well observed inside the detector. The solid red line is for 
the oscillation parameter values in the (3+2) model as given in the upper row of Table [2] 
The blue dashed line shows the results for the (3+1) best fit values: Am|^ = 1.78 eV^ and 
sin^2^ee = 0.089 obtained from reactor anti-neutrino data with the new predictions for the 
reactor flux 26 . As seen in Figjsj the shapes are quite different for (3+1) and (3+2) waves 



Also, a comparison between the amplitudes of the wave in various LjE bins cancels flux 
normalization and background systematic uncertainties to a large extent. 

In Table [8] we present the amount of v^. flux that is needed in 5 to 50 kt LENA type detectors 
to exclude the two different sets of oscillation parameter values in the (3+2) scheme as given 
in Table [21 The results are presented at 3 cr CL (4 dof) under the assumption that there are 



no beam or non-beam backgrounds for LENA disappearance search (see Sec. 4.3). In this 
analysis, we consider 15% systematic error on the signal and we follow the same numerical 
method used for the appearance results. Table [8] shows that the required amount of flux is 
highly dependent on the choice of the parameter values rather than the detector size. For 
example, we need ~ 26 times more flux for the parameter set A compared to B in the 25 kt 
LENA option. In Table |9| we present the results for NOz/A considering a 20% systematic 
error on signal and a 5% systematic error on Michel decay backgrounds as measured during 
beam-off running. From this table, one can see the impact of the different choices of effective 
background on the required flux as compared to the no-background case. 
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Figure 5: This figure shows the ratio of events with and without sterile oscillation as a function of 
the reconstructed L/E in 50 kt LENA. Here we have used total 4 x 10^^ v^- The solid red line is for the 
parameter values in (3+2) model as given in the upper row of Tablc[2] The blue dashed line is computed for 



(3+1) best fit values: Am|]^ = 1.78 and sin^ 26'ee ~ 0.089 as obtained from 26 using reactor anti-neutrino 
data with new predictions for the reactor flux. 



Fiducial Mass 


Flux 


Flux 




(A : Ref. [261 ) 


(B : Ref. f25]) 


50 kt 


3.39 X lO^'^ 


0.093 X lO^'^ 


25 kt 


5.55 X 10^0 


0.214 X 10^0 


10 kt 


11.25 X 10^0 


0.569 X 10^0 


5 kt 


22.1 X 10^° 


1.202 X lO^o 



Table 8: It shows the amount of neutrino flux needed to exclude the two different sets 
of oscillation parameter values as given in Table |2] at 3 a CL with 4 degrees of freedom 
(Ax^ = 16.25) using 5 to 50 kt LENA in disappearance mode. 
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Flux 
(A : Ref. |26|) 


Flux 
(B : Ref. (25)) 


50000 


5.9 X lO"^^ 


1.325 X 10^ 


25000 


4.615 X 10^^ 


0.963 X 10^1 


10000 


3.408 X 10^1 


0.636 X 10^1 





1.742 X 10^1 


0.0945 X 10^1 



Table 9: It depicts the amount of neutrino flux needed to exclude the two different 
sets of oscillation parameter values as given in Table [2] at 3 cr CL with 4 degrees of freedom 
(Ax^ = 16.25) using 14 kt NOz^A far detector in disappearance mode. We show the required 
amount of flux for four different choices of effective Michel electron background. 
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Figure 6: Expected constraints from DAR-LENA setup on sterile neutrinos in the (3+1) model at Scr 
CL (2 dof) using disappearance search. The dashed green curve shows the 99% CL (2 dof) limit from 



reactor anti-neutrino data with new reactor fluxes 66 . The triangle and the bullet show the (3+1) best- fit 



values of all reactor data with old and new fluxes respectively. Results are shown for 5 to 50 kt LENA type 
detector with 4 x 10^^ Vp. 
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Figure 7: Expected sensitivity limit of DAR-NOi^A setup to (3+1) sterile neutrino oscillation at Sct 
CL (2 dof) in disappearance mode. The dashed green line shows the 99% CL (2 dof) limit from reactor 
anti-neutrino data with new reactor fluxes 66 . The triangle and the bullet show the (3+1) best-fit values 
of all reactor data with old and new fluxes respectively. Results are presented considering 25000 and 50000 
effective Michel electron backgrounds with 4 x 10^^ (100 kW) and 4 x 10^^ (1 MW) Ve- 



in Fig.[6]we show the sensitivity hmit of the DAR-LENA setup to sterile neutrinos in a (3+1) 
model at 3 (T CL (2 dof) using the disappearance channel. We compare our results with the 
99% CL (2 dof) limit from reactor anti-neutrino data with the new reactor fluxes [66]. The 
triangle (Am^^ = 1.78 eV^ and sin^ 26'ee = 0.058) and the bullet (Amf^ = 1.78 eV^ and 
sin^2^ee = 0.089) show the (3-1-1) best-fit values of all reactor data with the old and new 
fluxes respectively. These are shown as benchmark points to judge the performance of the 
LENA setup. A 10 kt LENA with a flux of 4 x 10^^ Pe is sufficient to cover these test points 
and can provide stringent test of the recent reactor anomaly with high significance. Fig. [7] 
shows the sensitivity limit of the DAR-NOz/A setup to sterile neutrinos in a (3-1-1) model. 
Again, results are presented at 3 a CL (2 dof) for the disappearance mode. We show the 
results considering 25000 and 50000 effective Michel electron backgrounds with 4 x 10^^ (100 
kW) and 4 x 10^^ (1 MW) Ve total fluxes. For NOz/A, a 100 kW machine is marginal in 
covering the test points and a higher-power, full DAE5ALUS type machine, is needed to 
cover the entire parameter space. By comparing Fig. [6] and Fig. [7| one can say that a 14 
kt NOz/A with 10 times larger flux can give similar performance as to a 50 kt LENA with 
a reference flux of 4 x 10^^ Vf.- 



6 Summary and Conclusions 

A host of recent SBL neutrino oscillation experiments have provided hints that there may 
be oscillations to sterile neutrinos with a mass squared difference of the order 0.1-10 eV^. 
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Sterile neutrino models involving three active and one (3+1) or two (3+2) sterile neutrino 
states have been proposed to explain these results. These models demand that there be 
both appearance and disappearance associated through or to sterile neutrinos. Thus, to 
address the validity of these sterile neutrino models, better precision measurements of both 
appearance and disappearance oscillations in this mass region will be needed. 

The combination of a long liquid-scintillator neutrino detector combined with a cyclotron 
DAR neutrino source is very sensitive to neutrino oscillations for both the z/g appearance 
and the v,, disappearance channels. Such an experiment can observe the L/E variation 
of the oscillation rate within the detector and, therefore, provide proof that the data is 
explained by neutrino oscillations rather than other types of models. In addition, using the 
L/E variations within the experiment makes these types of measurements fairly insensitive 
to normalization uncertainties and backgrounds that do not have the L/E dependence of 
neutrino oscillations. 

The future LENA experiment provides an example of the capabilities of a very-large, liquid- 
scintillation detector. As shown above, the LENA detector provides unmatched sensitivity 
in the Ve appearance channel and completely covers the LSND and MiniBooNE signal regions 
at more than 5 a confidence level with a cyclotron power of 10 kW combined with a 5 kt 
detector. This would be a definitive investigation of the the LSND/MiniBooNE reported 
signal and the distinct L/E dependence of the appearance signal would provide unique 
confirmation that the signal is associated with high-Am^ oscillations. 

In addition, the LENA experiment could provide a robust search for z/g disappearance by 
observing the L/E dependence of the detected electron neutrinos. The sensitivity of the 
experiment would cover the region suggested by the recent reactor Pg disappearance obser- 
vation at 3 cr and, again, would provide key information on the possible interpretation of 
the result through the L/E dependence of the disappearance rate. 

On an earlier timescale, the NOi^A detector, now under construction, combined with a 100 
kW neutrino source could provide an initial test of disappearance at high Am^ with 
mixing angle sensitivity in the range of the recent reactor result. Since such a measurement 
is dominated by statistical uncertainties, making a NOi/A measurement with a 1 MW source 
would have a 3 cr sensitivity at small mixing angle down to sin^ 26ee ~ 0.02. 

In conclusion, we have shown that large neutrino detectors using liquid scintillator combined 
with high intensity 10-100 kW cyclotron DAR neutrino sources would have unprecedented 
sensitivity to sterile neutrino oscillations in the 0.5-10 eV^ mass squared difference region. 
Such experiments are therefore an important option as a next step in investigating neutrino 
oscillations to sterile neutrinos. 
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